PUPT 1480, lASSNS 94/51 
hep-th/9408XXXX 
August 1994 



Geometric Entropy, Wave Functional, and Fermions 



-k 

Finn Larsen 



Department of Physics 
Joseph Henry Laboratories 
Princeton University 
Princeton, N.J. 08544 



Frank Wilczek^ 



School of Natural Sciences 
Institute for Advanced Study 
Olden Lane 
Princeton, N.J. 08540 



■k Research supported in part by Danish National Science Foundation Fellowship. 

larsen@puhep 1 . princeton . edu 
t Research supported in part by DOE grant DE-FG02-90ER40542. wilczek@iassns.bitnet 



ABSTRACT 



We develop techniques for calculating the ground state wave functional and 
the geometric entropy for some simple field theories. Special attention is devoted 
to fermions, which present special technical difficulties in this regard. Explicit cal- 
culations arc carried through for free mass bosons and fermions in two dimensions, 
using an adaptation of Unruh's technique to treat black hole radiance. 
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1. Introduction 



Recently there has been great interest in the concept of geometric entropy [1-9] . 
Almost all of the explicit work has been done for free scalar bosons, although one 
can adapt conformal field theory arguments to cover some additional special cases. 
Among several possible generalizations, perhaps the most interesting is to fermions. 
Fermions play a role in several of the speculations regarding geometric entropy, 
especially those involving its analogues in superstring theory. The ground state 
wave functional and the path integral, fundamental objects from any perspective 
and certainly central to any discussion of geometric entropy, have quite a different 
character for fermions than for bosons. Whereas for bosons the states of the 
quantum field theory can be labelled by classical field configurations, since the 
field operators on a complete hypersurface constitute a complete set of commuting 
observables, the corresponding structure for fermions is less transparent. 

Let us recall the basic construction used to define geometric entropy. One 
divides space into two regions, which we shall call inside and outside. We consider 
a quantum field theory defined in all of space, and possessing a complete set of local 
observables ^{x). We divide these into two disjoint subsets according to whether 
their argument x lies on the inside or on the outside. Any state may be labelled 
by its amplitude as a function of the the eigenvalues (^in, i^out of the ^{x) with x 
lying on the inside or the outside respectively. Now consider a definite pure state 
described by the wave functional ^'[Ciiu Cout]- In the body of this paper, we shall 
be concerned with the ground state exclusively, and when we speak of geometric 
entropy without further specification we shall have this case in mind. Then the 
density matrix appropriate to an observer who only has access to the outside region 
is obtained by tracing over the variables localized in the in-region: 



This density matrix in general no longer describes a pure state, and one defines 
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the geometric entropy relative to the given state and the given division of space as 

5geom = -TVplnp . (1.2) 

This 5" is a measure of the strength of the correlations between the inside region 
and the outside region, information concerning which is lost in the process of trac- 
ing over the inside region. It has several qualitative properties that are rather 
different from those of thermodynamic entropy; notably it is not extensive, and its 
value is invariant under interchange of "inside" with "outside". In field theories 
with sufficiently singular ultraviolet behavior (that is, in essentially all ordinary - 
non-topological - field theories) the geometric entropy is dominated by short-range 
correlations. In the ultraviolet limit all theories look scale invariant, and further- 
more the important correlations arise from field fluctuations in the direction of the 
normal at the boundary. As a matter of technical convenience, rather than concep- 
tual necessity we shall, below, often use free massless fields in 1-1-1 dimensions as 
our working material. Due to the considerations just mentioned, this specialization 
is less severe than might appear at first sight. (In this regard, the interested reader 
might like to refer to Equation 36 of [4] as an illustration of the special significance 
of the 1-1-1 dimensional massless theories. Such theories, of course, also occupy 
center stage in superstring theory.) 

Although we will not pursue it much further in this paper, let us briefly note 
that one can also define the geometric entropy for mixed states, e.g. thermal states. 
In that context several of the points made above require modification. In particular, 
there is no longer symmetry between inside and outside; futhermore, the difference 
between the thermal geometric entropies at different temperatures is not expected 
to be dominated by ultraviolet behavior, and it should be approximately extensive 
for large volumes - the bulk geometric entropy density coincides with the bulk 
thermodynamic entropy. 

We have considered the geometric entropy in general conformal field theories 
previously, and obtained quantitative results which coincide with the ones pre- 
sented here [8] . The goal of the present work is, in a sense, to do it again the hard 
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way - working with field variables (especially, fermion field variables ) directly, 
without invoking conformal symmetry. Our goal is, by being very pedestrian and 
explicit, to sharpen our understanding of the technical issues involved, particularly 
in defining the wave functionals and the path integrals. 

As an example of the sort of subtlety that arises, consider the following. In 
their standard formulation, relativistic fermions obey a linear equation, arising from 
a quadratic action. One is accustomed to evaluating quadratic path integrals by 
substituting the solutions of the classical equations of motion (with the appropriate 
boundary conditions) back into the action. But here the action will simply vanish, 
and by proceeding naively one obtains a trivial - wrong - answer. 

Below, by exercising more care, we shall be able to obtain meaningful, useful 
expressions in this and similar contexts. Having these in hand, we will proceed 
to evaluate the density matrix and geometric entropy for free massless bosons 
and fermions in 1+1 dimensions. When the dust settles, we are left with a very 
convenient formalism that treats bosons and fermions on an equal footing. Here, 
as in most other work on geometric entropy, the insight that allows us to pursue 
matters to the very end is an inspired choice of variables, first used by Unruh [10], 
who used it in a different but related context. But there are some small surprises 
in the details: our generalization of the Unruh ansatz contains a funny coefficient, 
signs must be treated very carefully, and boundary conditions are nontrivial. 

2. The Wave Functional 

To discuss geometric entropy we need to describe our reference state in terms 
of local variables. In ordinary quantum mechanics we can project on eigenstates 
\q) of the position operator q, thus defining '${q,t) = (g|\I'(t)). For field theory, we 
want to choose a complete set of local variables generalizing q. The field operators 
0(a;) form a natural choice, but it is by no means unique. We might, for example, 
choose instead to diagonalize the canonical momenta 7r^(x). It is by no means 
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obvious at first sight what constitutes the "natural" choice is for fermions; we shall 
see it is in a sense a linear combination of the two just mentioned. 



2.1. SCHRODINGER EQUATION AND PATH INTEGRAL 

Although the definition of geometric entropy docs not refer explicitly to dynam- 
ics, we will find it useful to consider solutions to the time-dependent Schrodinger 
equation 

(2.1) 

where H is the Hamiltonian operator. The formal solution is 

\^^{tf)) ^Te-'^^*f-*'^\^{ti)) (2.2) 



where T denotes time-ordering. Time-ordering is an essential requirement, because 
the Hamiltonian is an integral over objects that do not mutually commute. 

The formal expression (2.2) for the wave function is difficult to evaluate exphc- 
itly for specific field theories, and the time-ordering obscures the Lorentz invariance 
of the theory. Indeed, the Hamiltonian as well as the time-ordering procedure de- 
pend explicitly on the choice of reference frame. Fortunately, one has the identity 

^[0/] = {(j)f\Te-'^^*f-*'^\(i)i) = j^^'P e*"^^^^ • (2.3) 

In this path integral, one integrates over all field configurations for which the 
field is a prescribed value 0j at some early time ti and equal to 0/ at the final 
time tf. The field configuration at the early time selects a specific state. The 
measure will be constructed below. If the action A((^) and the measure Vcj) are 
Lorentz invariant, the transformation law for the wave functional is simply given 
by appropriate change in boundary conditions. 



6 



2.2. Symmetry and the Geometric Derivation 



We recall a recent simple derivation [11] of (2.3), that is very instructive for 
our purposes. The Hamiltonian expression satisfies a simple first order differential 
equation. Since both expressions reduce to 1 for tf = ti, to prove their identity 
it suffices to show that the Lagrangian expression satisfies the same differential 
equation. Thus we must evaluate 



5 Jvcf) e 

where the variation is a change in tf to tf + 5tf. 

The definition of the energy momentum tensor T^,/ is 

dA^-^J T^.dg'^'d'^x (2.4) 

under a change dg^'^ in metric. The variation 5 may be implemented by choosing 
a new ^oo = (1 + '^'^~f + const) . Following [11] (using, in essence, a special case 
of the Schwinger action principle) , we require the expression to be invariant under 
changes in coordinates, thus deriving 

5 jv(t)e'^^^ j V(j)e'^ j T^.^dg^'^d'^x (2.5) 

for true changes in geometry. Notice the sign of this expression. Next, we can use 
conservation of (integrated!) energy-momentum to put this in the form 



SjV(j)e''^ _ I 1 



Stf 2ntf — tq 



J VcPi J TQQd'^x) e'^ = -m[ct)] , (2.6) 



which establishes (2.3). The action of Hamiltonian operator H on the wave func- 
tional is defined by the last equation, i.e. it is found by evaluating the operation 
occurring in (2.6) on the final shce of the functional integral. 
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For the preceding derivation, the crucial property of the path integral mea- 
sure employed is invariance under time dilation. We will also require invariance 
under translation with a function. A measure satisfying these properties will be 
constructed in the next sections, and will involve some surprising subtleties. 

To motivate the detailed investigation, let us elaborate further upon an dif- 
ficulty with with the naive application of the path integral for fermions, that we 
already touched on briefly. For simplicity, we consider a Weyl fermion in two di- 
mensions. On shell this can be represented as a function of (say) x + t, or in 
complex variables as a holomorphic function of z. The action is A = ^ J i/jBiIj. 
Now it appears that, at the saddle point, this action vanishes, since the equation 
of motion is simply dip — 0. One the other hand, consider the Hamiltonian expres- 
sion for the wave functional: the Hamiltonian H — J ipd^ip certainly does not 
vanish for holomorphic ■^s. How, then, can the two expressions (2.3) for the wave 
functional possibly agree? 

To address this question recall that in the derivation of the path integral for- 
mula for the wave function we exploited the feature that the measure is invariant 
under a rescaling of the time variable. But under such a rescaling, z mixes with z. 
It is therefore important that, as intermediate configurations in the path integral, 
we use a set of states that is closed under time translations. That means that 
even for chiral particles such as Weyl fermions, one cannot impose the equation 
of motion on intermediate states. To get a simple path integral we must allow 
both right and left movers in intermediate states, even for a Weyl particle. Once 
we do that, it becomes less mysterious why the action does not necessarily vanish. 
Note that in this reasoning we tacitly assumed that a stationary phase method is 
applicable and exact for the quadratic action, or in other words that by imposing 
the equations of motions found by varying the action with respect to the field, we 
find a field ipd that saturates the functional integral. In the next sections we will 
construct the path integral that has this desirable feature explicitly, and verify that 
indeed such a classical field exists even for chiral fermions - but the classical field 
will include both left and right movers. 
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2.3. The Holomorphic Path Integral 

A popular modern way to construct path integrals uses the holomorphic repre- 
sentation, also called the Bargmann-Fock representation [12-13]. In this represen- 
tation we work in bases that diagonalize the positive and negative frequency parts 
of the field, rather than the Hermitean operators and tt oc dt4> which are perhaps 
more natural from an intuitive standpoint. In the holomorphic representation we 
diagonalize the (f){x) + in{x) and thus also their Fourier transforms, denoted 0;^. It 
is in terms of the these latter operators that the free Hamiltonian can be expressed 
as a sum of fundamental harmonic oscillators, i.e. H = J In principle 

all we need for the calculation of the geometric entropy (or anything else) is the 
ability to evaluate a complete set of commuting local variables, so either choice 
will serve. 

For a boson field (p we use the coherent ket states 

10) = Yl e-^^^'^'^e^l'^'lO) . (2.7) 
k 

This state is an eigenfunction of the positive frequency part of the field operator, 
i.e. ^k\4>) — 4>k\4'k)- The index k labels the degrees of freedom which, for a 
free scalar field, can be identified with the momentum. For massless fields in two 
dimensions k may be restricted to take either only positive or only negative values, 
if one of these restrictions is imposed the field is chiral. The bras are obtained by 
complex conjugation of the kets. They diagonalize the negative frequency part of 
the field to the acting to the left: {4>\(pl. = {4>\4>k- There is a completeness relation 

I lld^kd<l>k = (2.8) 

k 

which implicitly normalizes the measure. 
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Now we can evaluate the evolution operator, as follows: 



n—l n—1 



k j=l j=0 

n—l n—l 

k j=i j=o 

Hill d^jd<pj) n e-f^(^^+^^^-^^^^-)+'"'=^^+^'^^-^]^*^- . 

k 3=1 3=0 

(2.9) 

In this manipulation, we handled the time ordering by inserting a sequence of 
complete sets and used the fact that coherent states remain coherent under time 
evolution. To avoid proliferation of indices, we suppressed indices k on (j) and 0; the 
remaining indices label the time slices, with the understanding that times and n 
replaces i and /. Taking Atj = tj+i — t j — > and using the notation = ^^Af""*^^ 
we find 

U{(f>f\(f>i) = J T>(j) [e5(^/'^/+'^i'^*)e~/'^^ [\{4>4>-H)+i'^4>'l>]^^Q-\{4>f4>f+<^>^<l>i) (2.10) 

in which, of course, we have introduced a new notation for the measure. The 
purpose of the additional, vacuous exponentials will be explained shortly. 

A Gaussian integral worthy of the name should be calculated exactly using 
the stationary phase method. Here a first cut at such an evaluation would involve 
imposing the equations of motion found by varying and separately*. How- 
ever to obtain a precise prescription the correct measure must be kept in mind. 
That measure included only variables with indices 1 through n — 1, so in apply- 
ing the stationary phase method we must not allow variations at the end points. 
In the given variables this is an awkward constraint, because as we vary with 



★ We ignore the determinant, which does not depend on the boundary conditions 
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(f)f fixed we inevitably vary 0j, and similarly for the holomorphic variables. By 
adjusting the boundary term appropriately, one is able to organize the expression 
as in (2.10), where the expression in the square brackets allows the full variation. 
Indeed, upon translating back to the discrete notation, one easily checks that the 
bracketed expression contains neither 0j nor 0^. Evidently the symbols with those 
names should be considered functions of the other boundary conditions, not as 
independent parts of the path integral. By way of contrast, in the factor outside 
the brackets, the bar really does denote complex conjugation. 

The expression in the square brackets is the standard expression for the evo- 
lution operator in the holomorphic representation. Finding and integrating the 
equations of motion, one readily calculates 

U{4>f\cl)i) = n [e^' --"'='*^-'^V.]e-i(^/0/+0.0O (2.11) 
k 

This expression could also be derived directly, using simple properties of coherent 
states. 

Taking the initial state to be vacuum, i.e. (f)i — 0, we find the vacuum 
wave functional to be a simple normalized Gaussian. This final result should 
come as no surprise, because the Hamiltonian is zero acting on the vacuum, so 
U{(/)f\0) = {(pflO), leaving only the normalization of the coherent state. This sat- 
isfying result exhibits also, unfortunately, the essential triviality of the procedure: 
the insertion of intermediate states, and subsequent integration, is vacuous. The 
entire calculation effectively occurs at one time slice, and one does not find the 
Lorentz or conformal symmetries exhibited explicitly. 
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2.4. Feynman's Path Integral 



To exhibit these symmetries, we would hke to find a form of the path integral 
that, written in real space, resembles the one derived from geometrical consider- 
ations. For this purpose, we return to the original Feynman construction of the 
path integral. The Hamiltonian is 

iJ = /g-4^»-/^^tfI + -M); (2.12) 

where 4 = ^M^Mj^ and ^J. = so that = 27r5(A; - k'). As a 

basis in the Hilbert space we introduce the overcomplete set |g) with the property 
Qk\Q) — From the commutation relations \pkAli[ — — 27ri5(A; — k') we derive 

that act on such states as ^fc|g) = i-^\q). Introducing also eigenstates of 
the momentum operator we find {p\q) = e~*^'=^'=. The operators % and are 
Hermitcan so their eigenvalues arc real. From the Hamiltonian we see that vacuum 
is defined by 0fe|O) = 0, and integrate to find 



= (g|0) cx JJe"^"^ (2.13) 



Normalizing the measure so J dqk e"'^^* = 1, 



J dq\q){q\^l^ J dp\p){p\ (2.14) 



(we omit the indices and denote the numbers q^ collectively as q, and similarly for 
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Pk), we are ready to derive the path integrah 



U{qf\qi) = {qf\ Te-'^qi) 




n—l n—1 



lldqjlldpj{qj+i\pj){pj\e '"^'^qj) 



j=l j=0 




n—l n—l 



(2.15) 



k 




In this expression, we have defined tn — tf, tQ — ti, and Atj — tj+i — tj, and used 



a continuum notation. 

We have gone into some detail here so that the main point comes out clearly and 
unambiguously, as follows. By construction, the variables q^ are real. This property 
must hold even it the field is chiral, that is if the k are restricted to be positive 
(or negative). The classical field includes both components nevertheless, because 
that is the only way it can be real. All this conforms with our earlier remarks on 
the geometric derivation of the path integral. We reahsed in that context that one 
must allow both left and right movers in the intermediate configurations, so as to 
work in a set that is closed under dilation of time. 

To evaluate the integral, we write 4> — 4>c[ + S4> where (pel is a solution of the 
equation of motion with the specified boundary conditions, while 5(f) vanishes on 
the boundary but is otherwise arbitrary. Inserting this decomposition in the path 
integral, one easily sees that the cross terms disappear. Hence the path integral 
factorizes, with only one factor depending on the boundary conditions. And that 
factor is simply the integrand in the path integral, i. e. exponential of the classical 
action, evaluated for the field <pci- 

It is easy to calculate the classical field even for general boundary conditions, 
i.e. (j) fixed to he (f)f at tf and (pi at U. We are especially interested in the special 
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case where the initial state is the vacuum. This case is most easily handled by 
noticing that taking tj — ti —ioo the Hamiltonian projects on to the ground 
state. Indeed, finding the classical field and taking this limit, with (f>i fixed at 
some finite value, has the same effect as removing the components that approach 
oo in the limit, or alternatively requiring that the classical field satisfies — > in 
this limit. Hence the wave functional for the vacuum state is calculated simply by 
imposing — > as t — > —ioo on the classical field. 

To tie up this discussion let us relate this careful explicit evaluation, with an 
identified measure Y\k ^Qk and a definite prescription for calculation of the path 
integral, to the heuristic (2.3). We shall illustrate this by reference to the wave 
functional for a chiral field in 1+1 dimensional Euclidean space. We find the 
classical field that satisfy the equations of motion and the appropriate boundary 
conditions, and calculate the path integral for that one configuration. The vacuum 
state is identified by the boundary condition — > as t — > —ioo or, in Euclidean 
space, ^ as T — > — oo. Imposing the reality condition, we write 

oo 

0ei = J ^== e'^' + ^ke'''] . (2.16) 

— oo 

Now we select the vacuum by restricting the integral to positive k. Calculating the 
action for this field configuration, we find 

oo 

If - 1 f dk - 
= ^ y 50ci^<^ci = 2 y ^40A: (2.17) 



The action integral was over the r < half plane. Finally we write the wave 
functional 

= ]^e-^^'=<^'= , (2.18) 
k 

a result that coincides with the one found above in other ways. 
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2.5. The Path Integral for Fermions 



With these experiences in mind, let us turn to the question of finding a useful 
prescription for the fermion path integral. For free fermions the Hamiltonian is 
again of the form H — ujo)a, but the creation and annihilation operators satisfy 
anticommutation relations {a, a^} = 1. Introducing hermitean variables q — ^^=^ 

and p — -^-^7^, we find {p, = 0. This relation does not allow for a realization 
wherein q is diagonal and p is expressed as a derivative operator. That would 
require a nontrivial (anti) commutation relation. Thus the road leading to the 
Feynman path integral appears to be closed, and we must fall back on a version 
of the holomorphic integral. Indeed, one can introduce a basis that diagonalize 
a and represent as an anticommuting derivative. This leads us to the holo- 
morphic representation of the path integral, with the additional subtlety that the 
holomorphic variables are anticommuting, so that in manipulating the expressions 
we must take care of their order. In the bosonic (second-order) case this form of 
the path integral was rather trivial and unsatisfying geometrically, but as we shall 
now demonstrate the situation is quite different in the fermionic (first-order) case. 

As usual, we consider a massless field in two dimensions. Introducing coherent 
states 

1^) = JJe-^^'='^'=e^^'^'=|0) (2.19) 
k 

we preserve the normalization of the state and the resolution of the identity 



/ 



IV'XV'I^I ■ (2.20) 

A; 
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Since the Hamiltonian remains H = J |||fc|'*Aj|;V'fe we find as before 

k j=l 

= /"p^ JJe/'^*^[^^-'^^l+^l^l'^V' (2.21) 
k 

= J Vi/j e~^^'^^'^'' V'(9t-aa)V _ 
In transforming back to real space we liave introduced 

+ (2-22) 

which by construction is real. We see that for fermions, unlike for bosons, the 
holomorphic path integral leads to the invariant form that we know from the ge- 
ometrical derivation. However here as before these expressions are a little subtle. 
In particular, in deriving and interpreting them one may be tempted to give up 
the reality condition on the field, in which case one must take boundary terms 
carefully into account in evaluating them, and a trivial integral is found. 

To evaluate the path integral without giving up the reality condition we proceed 
as for bosons and write the most general field ^l} — ip^l + ^"0) where ip^i satisfies the 
boundary conditions on while Sip is on the boundaries, that is the initial and 
final time slices. Finally, choose ipcl to be the unique real field that satisfies the 
boundary conditions, as well as the Klein-Gordon (not the Weyl-Dirac) equation. 
(Any field with the proper boundary conditions might have been used; the Klein- 
Gordon equation is imposed only to insure uniqueness!). Now note that in 

^ J iPidt - da)iP ^^J - aa)^ci + SiPidt - da)Si^ (2.23) 

mixed terms do not enter, due to the reality condition on the field and the boundary 
conditions on 5ip. The path integral measure Vip is invariant under translations 
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with a fixed function ip^l} so we might as well take it to be VSi/j. Then the path 
integral factorizes. The fiuctuation part is independent of the boundary conditions 
and can be omitted. We find that the path integral can be calculated by taking 
the action of the field configuration ifj^i, as advertized. At no point in the present 
derivation did we use a variational principle - specifically, we never found that we 
should impose the equation of motion on the classical field. 

We have found that the fermionic path integral is naturally expressed using 
real classical fermion fields to parametrize the initial and final states, thus putting 
them on an equal footing. The interpretation of such an amplitude between real 
states is not a priori clear, because the original problem did not have real fields. It 
is implicit in our derivation that the path integral expresses an amplitude between 
coherent states. To find the ket-coherent state that corresponds to a given real 
wave function, take the spatial Fourier transform. The components with positive 
wave vector provides the eigenvalues of the elementary oscillators of the Weyl 
field. Conversely, the real classical field corresponding to given eigenvalues of the 
elementary oscillators is found by using the eigenvalues as Fourier components with 
positive k- vectors, and their complex conjugates for the negative k- vectors. Thus 
the amplitudes that are expressed by our unorthodox path integral are exactly 
the same as those that are expressed by the more conventional holomorphic path 
integral. 

Just as for bosons, on taking ti — > —ioo with ipi fixed we find that il^i disappears 
from the problem and is replaced by the requirement that ipcl — > as t — > —ioo. 
This is a very convenient characterization of vacuum. As a final step in elucidating 
the fermionic path integral, let us Euclideanize the classical field and the action. 

To obtain expressions that resemble those found previously for bosons, we require 
that the Euclidean classical field is real; that is, that it is a sum of a holomorphic 
function and an antiholomorphic one. We then have, finally. 



It is easy to verify that the exponent is real. In this form, conformal invariance is 
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manifest. 

To exemplify the use of this machinery, let us calculate the vacuum wave func- 
tional. We proceed as for bosons: write the classical field 

oo 

^ei = y ^ [Vfee^^ + Vfee^^l (2.25) 



where we have imposed the reality condition, and the condition that the field 
vanishes at early times. Inserting this in the action we find 

If - 1 f dk - 
AiiJcl) = ^J V'cl^^cl = 2 y ^^fc^fc (2-26) 

and for the wave functional 

= n 6"^^'=^'= . (2.27) 

A;>0 

Thus this formahsm indeed yields the expected (O|'0)- 

3. Geometric Entropy of Free Bosons and Fermions in 2 Dimensions 
3.1. Strategy 

We now return to the problem that motivated the preceding ordeal. We will use 
the formalism developed to carry out the calculation of the geometric entropy for 
massless scalar bosons and spin ^ fermions in 1+1 dimensions. In this context, the 
essential problem is the entropy associated with a half-line. We will demonstrate 
that each step in the calculation of the entropy can be carried out very explicitly, 
using the flexibility of the path integral expressions just developed. This may not 
be the most efficient way to reach that specific goal, but it provides an explicit, 
and hopefully transparent, derivation of the entropy in a manner that is parallel 
for fermions and bosons, and capable of generalization. 
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The first step is to calculate the vacuum wave functional 

*[i?,L] <x Jv(f) e-^^"^) . (3.1) 

Here the path integral is over all fields that satisfy 

(f)(a) = e{a)R{a) + 9{-a)L{a) (3.2) 

at a time slice taken to be r = 0, and we also impose 

4){x) r — > — oo 

to project onto the vacuum. In earlier sections we worked with the Fourier com- 
ponents of the field but these are not locahzed to be either on the left hand side 
of the axis, or on the right hand side. Hence we need to choose as a basis instead 
functions that are partly localized, but still resemble Fourier modes sufficiently 
to approximately diagonalize the action. Wavelets [14] are designed for exactly 
this purpose, that is to provide wave-like functions with compact support. They 
diagonalize the action approximately, and are likely to be of considerable use in 
problems with more complicated structure*. For the present however, we will stick 
to Fourier modes, and simply transform the argument of the fields R and L in- 
stead. More precisely, we will introduce a convenient coordinate system that maps 
both half lines to full lines, for which we can use the standard Fourier transform. 
The density matrix will not quite diagonalize, but will break up into 2x2 blocks. 
This trick is in essence due to Unruh [10]. The boundary conditions have a unique 

★ There is also a more fundamental point that ought to be mentioned in this context. The 

geometric entropy as defined corresponds to the density matrix for a hypothetical experi- 
menter who has complete access to arbitrarily high frequency modes on the outside, but no 
access to the inside. A more realistic idealization would be to allow access to low-frequency 
modes on the outside only, tracing over both very high frequencies in general and also low 
frequencies on the inside. These notions could be formalized using wavelets. The entropy 
thus defined would be finite, diverging only as the limiting frequency is taken to infinity. 
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solution among fields that arc the sum of a holomorphic and an anti-holomorphic 
piece. That field is the classical field and the path integral is calculated by find- 
ing the action of the classical field. This is true, as we have seen, for bosons and 
fermions alike. 

Having obtained the wave functional in a convenient basis, the next step is 
to sum over the left variables, and then to find the entropy corresponding to the 
resulting density matrix. This will be done using a replica trick, as in [4]. 

We will use the complexification z — a + ir and z — a — ir. Since r is the 
Euclidean time this amounts simply to the light cone coordinates. This conven- 
tion interchanges r and a compared to the one conventional in the string theory 
literature. 

3.2. Classical Fields 

First we calculate the classical field for bosons. (f){x) is specified at r = 0; our 
task is to determine (f){x) in the entire lower half plane. We write 



oo 



(p{z,z) = 




— oo 



In the integral w is a real variable. This integral equation is clearly valid on the 
real line, and extends by regularity to the entire negative half-plane (it is just the 
usual Poisson integral for this problem) . We change of variables according to 

w — sign(w) e^, z — , z — , (3.4) 

leaving the field untouched. The field is defined in the lower half plane, so in 
inverting 2; = e'' we must choose Imr) < 0. Thus the positive half-axis is mapped to 
the entire real axis, and the negative half-axis is mapped to the line with imaginary 
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part —iir. We write 

R(x) = J ;^7=e"'""r, m = / ' (3.5) 

thus parametrizing the functions R and L by their Fourier components in the 
transformed variable. The reahty condition on the field is expressed as r^^ = f-^^ 
and = l-ui- Now we have 

oo 
— oo 

and calculating the integrals over x by contour integration and recalling imry < 
we find 

It is easy check that indeed (3.2) is satisfied, i.e. that (f) — R ior r] E TZ and 
(j) — L ior rj & TZ — in. Indeed, this expression could easily have been found by 
writing the general form of the wave function and determining the coefficients from 

the boundary conditions. In this reasoning, the boundary condition that — 
at early times, that is the choice of vacuum, is expressed by imposing regularity 
throughout the strip. The present, constructive approach has the advantage that 
it is easily generalized to the case of fermions. 

Indeed, let us write 

oo 

i;{z,z) = ^ f dw{— —)^{w) . (3.7) 

ZTT J W — Z W — Z 

— oo 

Let us again introduce the left/right split (3.2), and the change of variables (3.4). 
Now, however, we transform the fermion field according to 

i^{z,z)^e-'^''iP{r],f}), z^e'^,z^e^. (3.8) 

The necessity to transform V^, in contrast to (f), ultimately reflects the non-trivial 
conformal weight of ip. The original was real and that transforms to a real on 
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the real axis but to a purely imaginary on IZ — m. Introducing Fourier transforms 



this is expressed by f^^ = r_j^ and l,^ — —l-uj- We have chosen a different normal- 
ization here than in (3.5) for bosons, in order that the Fourier components r^^ and 
l,^ have mass dimension — ^ for fermions, as it did bosons. Collecting formulae, we 
have 

oo 

duJ f dx iojx —x[( ^^uJ 



iP(r],fi)e—^'^ / -;= / ^e-^'^^e^^ff— - (r] ^ fj)] 

f^yhu ' /27r J 27r ^^e^-e^ e^ + e*?^ ^' '^^ 



The extra factors of ea* and e^^ compared to the boson case come from the trans- 
formation of the fermion fields as a (i,0) field. This is also the origin of the extra 
i in the second term. Upon performing the integrals we find 

In principle this decomposition could have been found by writing the ansatz, 

and determining the coefficients from the boundary conditions. However this ansatz 

T) — fj 

is non-trivial, the e~2~ being due to the conformal dimensions of the fields. In 
the present, constructive approach it is well motivated from the transformation 
properties of the fermion field. 
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3.3. Wave Function als 

Having found the classical field for both bosons and fermions we proceed to 
find the wave function. The classical actions are 

oo 


and 

CX3 



In the fermion case wc avoid the notation to prevent confusion due to the relation 
lijj = —l-uj- In evaluating these expressions we have used the integrals 

The wave functionals are simply 

^[R,L] = e-^=' 

It is a good check, to verify some of their necessary qualitative properties. In the 
limit u ^ oo the wave functionals reduce to the results found with no left/right 
split, that is, Gaussians with the same normalization. This occurs because excita- 
tions that are almost localized do not mix with excitations on the other half line. 
At finite frequencies, there is overlap. We note however that the operator in the 
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exponent has the same determinant as previously, as indeed it must for a unitary 
change of basis. This feature is easily checked using the path integral measure 
Yi^ dl-^dlojdr-ajdrco, and the useful integral formula 

J dzdz e-''Mz+-zj+jz ^ (detM)^^ e^^~'^ 
where the upper and lower signs refer to bosons and fermions respectively. 
3.4. Geometric Entropy 

To calculate the entropy from the wave functional we must first calculate the 
density matrix 



p[R,R'] = y *[i?,L]*[i?',L] . 



The measure in the integral is simply Yl^ dl-ujdlcj, which we have normalized so 
that a Gaussian gives unity. The wave functional \E' is normalized by requiring 
Trp = 1 with respect to the same measure. We find 

-^t^' = n ^"^p^-2ihL;['^2^^^'"'^'' + '"^''^ " ^'^'^"^ + "^"'^^^^ ^^-^^^ 

for fermions as well as for bosons. We have defined Ir^^p = r-^r^j, and use again 
r_(^ = f,^. In verifying this expression it is important to recall that integrals 
over Grassmann variables gives a determinant in the numerator rather than in the 
denominator. 

Next we want to use the replica trick 

5'geom = -(1 - ^)lnTrp^ (3.11) 
so we need to calculate p^''. The result is 

^ n |^^^P^-^^hL^[^^2™(|r.|2 + \r'J) - (f.rL + f>.)]} 



(3.12) 
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where 



(2shn7ra;)2 



(bosons) 



(2chn7ra;)^ 



(fermions) 



(3.13) 



(2ch7ra;)2« 



These formulae can be verified inductively. The difference between bosons and 
fermions is two-fold. First: Grassmann integrals, as mentioned above, give deter- 
minants in the numerator rather than in the denominator. Second: in taking the 
trace for bosons we simply identify r,^ = and do the dr-^^dr^ integral, but for 
fermions we must take r^j = — instead. This difference has been explained in 
an elementary way by Soper [15]. We can also understand it simply in our frame- 
work, as follows. In the holomorphic formalism the typical bilinear operator can 
be expanded as a string of variables in the form anti-holomorphic, holomorphic, 
anti-holomorphic, etc. In taking the trace, however, we pair the last holomorphic 
variable in a string like this with the first anti-holomorphic variable. For Grass- 
mann variables this operation must be accompanied by a change of sign, which is 
most easily handled by changing the sign on one of the variables. This accounts 
for the antisymmetric boundary conditions. 

With these expressions the replica trick (3.11) can be carried through to yield 



= ±2/ ^(1 - T = 4 / d.^^^ (3.14) 



for the entropy. The upper sign refers to boson and the lower to fermions. These 
are simply the thermodynamic expressions for the entropy of a 1-dimensional gas 
of (spinless) bosons or fermions respectively. So far, the possible values of u have 
been specified informally as / and strictly speaking expressions such as Ylu; 
did not make sense without regularization. At this point wc impose periodic (or 
anti-periodic) boundary conditions in a box of length L, thus arriving at 



'geom 




(3.15) 



with c = 25 1 fermions and bosons, respectively. In this formula the length L is 
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the length as measured in the transformed coordinate system. Transforming back 
to the original coordinates and introducing an ultraviolet cutoff e and an infrared 
cutoff S we can write L = In — . This result and the interpretation of the ensuing 
divergence have been discussed extensively from another point of view in [8]. 

We have now accomplished the technical task we set ourselves, to obtain a path 
integral sufficiently ffexible to allow us to calculate wave functional and geometric 
entropy in a straightforward manner, applicable both to bosons and to fermions. 
We were surprised, that to do so we had to forge some new tools. 

Acknowledgements: We wish to thank Curt Callan, who was involved in the genesis 
of this work, for many helpful discussions. 
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